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Abstract We present a case study of the 13 July 2004 solar event, in which 
disturbances caused by eruption of a filament from an active region embraced 
a quarter of the visible solar surface. Remarkable are absorption phenomena 
observed in the SOHO/EIT 304 A channel; they were also visible in the EIT 
195 A channel, in the Ha line, and even in total radio flux records. Coronal and 
Moreton waves were also observed. Multi-spectral data allowed reconstructing an 
overall picture of the event. An explosive filament eruption and related impulsive 
flare produced a CME and blast shock, both of which decelerated and propagated 
independently. Coronal and Moreton waves were kincmatically close and both 
decelerated in accordance with an expected motion of the coronal blast shock. 
The CME did not resemble a classical three-component structure, probably, 
because some part of the ejected mass fell back onto the Sun. Quantitative 
evaluations from different observations provide close estimates of the falling 
mass, ~ 3 • lO""^^ g, which is close to the estimated mass of the CME. The falling 
material was responsible for the observed large-scale absorption phenomena, in 
particular, shallow widespread moving dimmings observed at 195 A. By contrast, 
deep quasi-stationary dimmings observed in this band near the eruption center 
were due to plasma density decrease in coronal structures. 



^ Institute of Solar-Terrestrial Physics SB RAS, Lermontov 
St. 126, Irkutsk 664033, Russia email: grechnev@iszf.irk.ru 
email: uralov@iszf.irl<.ru 

^ P.N. Lebedev Physical Institute, Lenisky Pr., 53, Moscow, 
119991, Russia email: slem@lebedev.ru 

^ Pushkov Institute of Terrestrial Magnetism, Ionosphere 

and Radio Wave Propagation (IZMIRAN), Troitsk, Moscow 

Region, 142190 Russia email: ichertok@izmiran.ru 

^ Ussuriysk Astrophysical Observatory, Solnechnaya St. 21, 

Primorsky Krai, Gornotaezhnoe 692533, Russia email: 

kuzmenko_irina@mail.ru 

^ Nobeyama Radio Observatory, Minamimaki, Minamisaku, 
Nagano 384-1305, Japan email: shibasaki@nro.nao.ac.jp 



(T].('ciiii(*\" ct. al. 



Keywords: Corona, Radio Emission; Coronal Mass Ejections, Low Coronal 
Signatures; Coronal Mass Ejections, Initiation and Propagation; Prominences, 
Active; Radio Bursts, Microwave (mm, cm); Surges; Waves, Shock 

1. Introduction 

Some major flares and coronal mass ejections (CMEs) are accompanied by 
large-scale disturbances moving over distances comparable with the solar radius 
like Moreton waves, weakly bright coronal transients (Thompson et al., 1998) 
usually termed 'EIT waves', and temporary depressions of soft X-ray (SXR) 
and extreme-ultraviolet (EUV) emissions (termed dimmings) appearing near 
eruption centers and far from them. They have been observed with the Extreme- 
ultraviolet Imaging Telescope (EIT; Delaboudinicrc et a/., 1995) on SOHO, the 
Transition Region and Coronal Explorer (TRACE; Handy et al., 1999), and other 
telescopes. Nevertheless, nature of all these phenomena is still controversial. 

Moreton waves are sometimes observed in the Ha line (Moreton, 1960) and, 
more recently, in the Hei 10830 A line {e.g., Vrsnak et al., 2002; Gilbert et 
al., 2004a). Uchida (1968) proposed them to be due to coronal disturbances. 
Balasubramaniam, Pevtsov, and Ncidig (2007) also argued their coronal nature. 
Moreton waves run faster (400-2000 kms~^) than 'EIT waves' (< 800 kms~\ 
typically ~ 250 kms"^). Warmuth et al. (2001, 2004a, 2004b, 2005) found their 
kinematical closeness and proposed that both chromospheric and coronal signa- 
tures of the waves could be created by a single decelerating disturbance, e.g., a 
blast shock (see also Hudson and Warmuth, 2004). Gilbert et al. (2004a, 2004b) 
found 'EIT waves' to be co-spatial with their counterparts in Hei. 

Some observations favor the wave nature of 'EIT waves' (Ballai, Erdelyi, and 
Pinter, 2005). A few moving bright features were interpreted as CME frontal 
magnetic structures (Derc et al., 1997; Uralov, Grechnev, and Hudson, 2005). 
Several authors propose explanations for 'EIT waves' within scenarios suggested 
for dimmings. Delannee and Aulanier (1999) proposed 'EIT waves' to be run- 
ning boundaries of successively opening magnetic structures ahead of expanding 
dimming areas, from which plasmas evacuate. Attrill et al. (2007) proposed to 
explain 'EIT waves' along with expanding shallow dimmings in their scenario 
based on interchange rcconncction. However, 'EIT waves' are not commonly 
associated with dimmings {e.g., Chertok and Grechnev, 2005b; Grechnev et al., 
2005). This means that their nature is not always the same. 

Trying to reconcile contradictory views, Zhukov and Auchcre (2004) proposed 
that 'EIT waves' represent both wave-like and eruptive components. Simulations 
of Chen et al. (2002, 2005) confirmed a possibility of two types of coronal wave 
phenomena during an eruption, one according to the scenario of Delannee and 
Aulanier (1999), and a fast coronal counterpart of the Moreton wave ahead. 

Quasi-stationary dimmings develop from ~10 minutes to several hours, reach 
—80% (drop up to 20% of the pre-cvent brightness; Chertok and Grechnev, 
2003a), and live up to two days (Hudson and Webb, 1997). Having once appeared 
in some region, they persist there (Chertok and Grechnev, 2005b; Grechnev 
et al., 2005). Such dimmings are related with active regions and large-scale 
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magnetic fields (Chertok and Grechnev, 2005b; Slemzin, Kuzin, and Bogachev, 
2005; Slemzin, Grechnev, and Kuzin, 2006; Zhang et ai, 2007). Their main 
interpretation is plasma density depletion due to its outflow from previously 
closed structures (see also Harra and Sterling, 2001; Harra et al, 2007). 

Some moving shallow darkenings seem to accompany 'EIT waves' in their 
expansion (see, e.g., Zhukov and Auchere, 2004). They can have the same wave- 
like nature as such 'EIT waves'. Another kind of short-lived, moving darkenings 
could be due to absorption of the background emission by material of ejecta. 
Such phenomena have been observed in images produced with EIT in the 195 A 
(Fexii) band and TRACE in the 195 and 173 A (Feix) bands (the 173 channel 
is usually called 171 to be consistent with the terminology of EIT; wo use 173 
to distinguish them). In particular, darkenings sometimes observed with EIT in 
the 304 A band (mainly Hen, 20000-80000 K) without counterparts in coronal 
bands (Chertok and Grechnev, 2003b) are expected to be certainly due to absorp- 
tion in the ejected filament material (Delaboudiniere, 2005). Such a phenomenon 
was recorded with the CORONAS-F/SPIRIT EUV telescope (Zhitnik et al, 
2002) in a major event of 18 November 2003 (Slemzin et al., 2004; Grechnev et 
al, 2005). A dark feature as large as the solar radius was observed at 304 A only 
to move during 1.5 hr across the Sun with the plane-of-sky speed of ~ 200 kms~^. 
Thus, the increasing observational material suggests that different kinds of 'EIT 
waves' and dimmings probably exist (as also Warmuth et al. (2001) proposed). 

'Negative bursts' (transient decreases of the total radio flux below the rel- 
atively quiet level of the solar emission) are believed to be mainly also due 
to absorption {e.g., Covington and Dodson, 1953; Maksimov and Nefedycv, 
1991), although other reasons are possible {e.g.. Sawyer, 1977). Since opacity 
of an absorber depends on frequency and plasma parameters, multi- frequency 
observations of such bursts promise quantitative diagnostics of ejecta. It seems 
to be useful to compare observations of an event, in which absorption occurs in 
both EUV and radio ranges. 

This paper addresses the 13 July 2004 event, in which a multi-component 
eruption, CME, Moreton and coronal waves, large-scale EUV absorption phe- 
nomena, long-lived dimmings, and a 'negative radio burst' were observed simul- 
taneously. The event included an impulsive M6.7 flare (00:09-00:23, all times 
hereafter are UT) in active region 10646 (N13 W46) and a CME with a central 
position angle of 294° observed after 00:54 with the Large- Anglo Spoctroniotric 
Coronagraph (LASCO; Brueckner et al, 1995) on SOHO. EIT observed large- 
scale disturbances in the 195 A band, in particular, a faint oval front at 00:24. 
A type II burst at 00:16-00:43 was reported by a few observatories. At 01:19, 
a large-scale darkening was observed with EIT at 304 A. During this event, a 
strong impulsive radio burst was recorded at 1-80 GHz followed by a decrease 
of the radio emission. The 'negative burst' was revealed in total flux records of 
the Ussuriysk Astrophysical Observatory (UAFO, 2.8 GHz), Nobeyama Radio 
Polarimeters (NoRP; Nakajima et al, 1985; Torii et al, 1979), and the Lear- 
month Observatory. The event was also recorded in the Ha central line in the 
Big Bear (BBSO) and Mauna Loa (MLSO) Solar Observatories, in the Hei line 
in MLSO, and in the TRACE 173 A channel. The Nobeyama Radiohehograph 
(NoRH; Nakajima et al, 1994) observed some fragments of the ejecta at 17 GHz. 
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Using multi-spectral data, we study this event and reconstruct a picture of the 
eruption. The ejecta disintegrated. One part fell back onto the Sun, while another 
part escaped. Dispersed material was responsible for the absorption phenomena 
observed. Estimations of masses of both these parts, performed using different 
methods, quantitatively confirm this picture. 

The event started with an explosive filament eruption and impulsive flare that 
likely produced a blast wave. Kinematics of the wave and ejected fragments are 
compared with analytical expressions, which match observations. Other aspects 
of the event concerning radio bursts, pulsations, and their association with the 
flare, ejecta, and shocks are addressed by Pohjolainen, Hori, and Sakurai (2008). 

Section 2 presents observational data. All measurements are related to the 
plane of sky (POS), if otherwise not specified. In Section 3, we address absorption 
in different emissions and estimate the mass of the ejecta. In Section 4, we discuss 
kinematical characteristics of the observed phenomena and compare them with 
theoretical expectations. In Section 5, we discuss results of our analysis, the 
overall picture of the event, and nature of 'EIT waves' and dimmings in this 
event. Section 5.4 summarizes our findings and their implications. 

2. Observations 

2.1. Filament Eruption, Moreton Wave, and Surge in Ha Images 

The 13 July 2004 event was observed in the Ha line center in BBSO and MLSO. 
No Ha flare was reported for this event, but images from both observatories show 
it and related eruptive activities. Figure 1 presents Ha images obtained in BBSO. 
The electronic version of our paper contains a movie 2004-07-13J3BSO_Ha.mpg, 
whose left panel shows Ha images with a limb darkening removed, and the right 
panel shows the same images after additional processing to reveal features of 
interest. Before the eruption, the active region contained a system of rather 
small filaments. Filament located between positions 1 and 2 starts to rise at 
00:02:30 with the top part accelerating at ~ 16ms~^, which reaches 7 Mm at 
00:13:25 with a speed of > 22kms-i (2"^ order fit). It erupts between 00:13:25 
and 00:15:25. A Moreton wave runs from the eruption site north/northeast (see 
Figure 2). At about its onset, surge 3 starts to grow from the western end 1 of the 
filament with the average speed of c:i 270 kms'^ (00:18-00:29). At 00:17-00:18, 
a surge from its end 2 appears and takes a V-like shape by 00:26 (frame c). 

The western surge 1-3 consisting of a multitude of thin, likely twisted, threads 
intensively develops. After 00:22, it starts to bifurcate into the northwest (1-4) 
and north (1-3) branches. Branch 1-4 rises up to its maximum POS height of 
~ 90 Mm by 00:30 and then mainly disappears. Branch 1 3 rises and slightly 
displaces east. New threads appear. Along with fragments emanating from the 
eastern vicinity 2 of the main filament and a new remote feature 5, they form 
to 00:30 a large-scale filamentary structure. Its components occupy a significant 
area 1-4-3-6-5-2 in frame d. Then this area continues to expand north and east, 
especially, due to the northern part of branch 3 and fragment 6 (frame e). Feature 
6 moves toward the North Pole with a speed ~ 155 kms~^ (00:31-00:53). As 
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MLSO Hei 10830 A images show, the disturbance probably associated with this 
moving feature reaches vicinities of the polar coronal hole at ~ 01:20. At about 
00:40, portion 3a of surge 1-3 overturns at a POS height of 190 Mm and starts 
to fall, while its farthest part 3 keeps on moving to the limb. After 00:50, the 
surge partially falls back and becomes fainter. This is due to either spreading out 
its components that decreases its optical depth, or their increasing line-of-sight 
velocities that results in the Doppler brightening. Throughout development of 
the surge, many threads and fragments keep connection with the position of the 
initial filament(s) 1-2, and their bases rest on its (their) vicinities. 

Prom Ha data we conclude that a filament eruption resulted in formation of 
some large-scale 'dome' of filament's threads and fragments to cover in several 
tens of minutes a significant northwest part of the Sun. 

2.2. High-Resolution TRACE 173 A Images of the Eruption 

TRACE observed the event starting from 00:07:27 in the 173 A channel mainly 
with intervals from 20 to 70 s, except for the gap of 00:21-00:48. Figure 3 presents 
several images of the eruption region, which TRACE permits to discern like a 
magnifier (see also movie 2004-07-13_TRACE_173.mpg). Pre-eruption images 
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Figure 2. Morcton wave (the front indicated by the arrows) in BBSO Ha images with 
subtracted average image. The dark wide ring is an artifact (a 'trace' of the active region 
that appeared in a toroidal Fourier transform in the course of background subtraction). All 
images are displayed nonlinearly to emphasize the wave. The cross marks the eruption center. 



show a dark filament (s) and a bright rope along it as well as the northern and 
southern loop systems surrounding the region. 

Four enlarged images (this area is outlined in panels e and / by dotted frames) 
in panels a d show two filaments located side-by-side to 'merge' and then to 
expand. The expansion of the combined filament is visible in frames e—l. The 
rising filament gets semi-transparent, and its parts become detectable {j,k). Loop 
LI expands in frames i-k and disappears by 00:18 (l). A flare starts beneath 
the filaments by 00:13:37. In frame i, bright, apparently diffuse material appears 
above the whole combined filament as its upper envelope, and later, in frames j- 
l, its brightest part appears as a jet-like feature Fl moving north (up). Another, 
the slower bright jet-like feature F2 moves west-northwest in frames j-l. Crossly 
dark/bright patterns are due to interference on the CCD matrix. 

The POS speed of feature Fl, v = 400 — 500 kms^^, is acquired in < 2 min 
(acceleration > 8 kms~^). With this speed and exposure times t^xp specified in 
frames, apparent extents of this feature are vte^pRQ'^^'^'^ / Rq^ — 18" — 31". Its 
jet-like appearance is caused by blurring due to its displacement dtiring the 
exposure. The same reason explains the semi-transparent appearance of the 
combined filament. It remains dark in all frames, where it is visible (by 00:17:48). 
The speed of its top is ~ (0.5 — 0.7)v. 

Hence, the eruption occurs between 00:14:28 and 00:15:20. The eruptive fila- 
ment remains cool. It is surrounded by a hotter envelope with probable temper- 
atures of 0.3-1.6 MK (the temperature sensitivity range of the 173 A channel). 
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Figure 3. TRACE 173 A images of the eruption: a— I, the early stage; m-s, the late stage. 
Insets a-d show enlarged images of filaments outlined in panels e and / by dotted frames. 
Labels LI in panels e—l mark an initial position of an erupting loop. Labels Fl and F2 
mark bright jet-likc moving features (position of label F2 is fixed). Insets p-s show enlarged 
areas outlined in panels m and n by dashed frames, in which a flow of absorbing material is 
indicated by the arrows. Another flowing absorbing feature is outlined by the ellipse in panels 
m and n. Observation times and exposure durations are specified in frames e— n. All images 
are nonlinearly displayed. 
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Expansion of the filament and adjacent structures starts before a motion of 
loop LI. Its displacement visible in frames when the ejecta was far from the 
loop, indicates passage of a wave. Thus, loop LI is pushed by the wave rather 
than vice versa. The first manifestation of the wave in frame i is somewhat ahead 
of the Moreton wave (Figure 2a). The speed of the loop is ~ 350 kms~^. 

Images obtained after the gap in observations show a post-eruptive arcade. 
Interesting are dark moving flows occulting brighter structures behind them. 
Examples are outlined by dashed frames in panels m,n (insets p-s show en- 
larged images, where a flow is indicated by the arrows) and by ovals. These 
flows represent moving absorbing (probably cool) material directed likely along 
magnetic field lines. Note also shrinkage of loops that we do not discuss here. 

The eruption destroyed the northern loop system. Loop LI erupted, but not 
reappeared afterwards. The southern loop system remained almost unchanged. 
The eruption looks like a 'directional explosion' pointed north. 

2.3. Coronal Disturbances in SOHO/EIT 195 A Images 

EIT produced 195 A images every 12 min (CME Watch program). Figure 4 
shows a pre-event (a) image and a late-stage (e) one along with six fixed-base 
differences h-d and /- h with the solar rotation preliminarily compensated (see 
Chertok and Grechnev, 2005b). Figure 5 provides additional information about 
large-scale faint disturbances, which are easier to see in running difl^erences. 

The eruption causes significant disturbances in the corona during 00:24-01:26, 
both brightening and darkening, that cover a huge area exceeding the NW quad- 
rant of the solar disk. A nearly perfect oval faint brightening is seen at 00:24 
(Figures 4c and 5b). It is not uniform in brightness, mainly 5-15 counts/pixel, 
while the brightness in quiet Sun's regions is about 40-50 and 5-21 in coronal 
holes. Three wide brightening sectors directed east-southeast, northeast, and 
west-southwest spread toward the oval from inside. The oval expands; its next 
position on the solar disk is denoted by a dashed arc in Figure 5c (00:36). 

Darkening areas are pronounced in Figure 4 (c,d, f-h). Some of them, mostly 
short-lived once, coincide with dark features like surges visible in Ha images {1- 
6, the same numbering). Some others are different, e.g., 7 and 8 and, especially, 
deep, long-lived dimmings. They show up as two regions Dl and D2 near the 
post-eruptive arcade, a dark arch between their tops, and a V-like dark feature 
on top of the arch. Comparison of panels h, a, and e shows that these deep, 
long-lived dimmings are due to significant darkening or disappearance of rather 
compact coronal structures, which were previously bright. Dimmings of shorter 
lifetimes with counterparts in Ha images are probably due to absorption of the 
195 A emission in the ejected material. A shallow depression (SD) moving toward 
the North Pole is faintly visible in fixed-base difference images (Figure 4d,g,h). 

The field of view of TRACE is insufficient to cover dimmings Dl and D2. 
TRACE 173 A images in Figure 6 show that the arch- like dimming is due to 
eruption of loop LL Dimming D2 is due to changes of the loop system near 
the western leg of the filament. The TRACE 173 A and EIT 195 A bands are 
different that explains some dissimilarity of their images. 
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Figure 4. SOHO/EIT 195 A images of tlie event, (a) Pre-event image, (e) later-stage image. 
Fixed-base (00:40 UT) differences (b— d and f-h) show large-scale coronal disturbances. Non- 
linear display is applied to emphasize them. Dashed black-white circles denotes the solar limb. 
Dashed white oval traces the aureole discussed in the text. Labels 1-8 denote darkenings; Dl 
and D2, long-lived dimmings. 
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Figure 5. Large-scale disturbances in the NW quadrant visible in running-difference 
SOHO/EIT 195 A images (nonlinearly displayed). The cross marks the eruption center. 
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Figure 6. TRACE 173 A images before (a), after (c) the eruption and their difference (b). 



2.4. Large-Scale Darkening in the SOHO/EIT 304 A Channel 

SOHO/EIT observed a large-scale darkening at 304 A in a single frame at 01:20. 
Figure 7a-c shows three sequential full-disk images produced in the 304 A 
channel. They are separated by 6 hours that rules out a possibility to study 
development and motion of the dark feature. It is certainly absent in the first 
(a) and last (c) frames, being well pronounced in frame b. Figure 7d shows 
a difference of frames b and a. Both positive and negative pixel values are 
restricted. The area bounded by the —25% level is 6.7% of the solar disk, and the 
deepest depression reaches —60%. The configuration of the large-scale darkening 
at 304 A is completely different from those observed in the 195 A band. 

2.5. Ejecta Observed with NoRH at 17 GHz 

The flare produced a strong microwave burst at 17 GHz, with its main phase 
lasting until 00:19 and a maximum brightness temperature up to Tb = 1.2-10^ K. 
To reveal faint moving features, we produced enhanced-sensitivity NoRH images 
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Figure 7. The northwest quadrant of EIT 304 A images produced before (a), during (6), 
and after the event (c); (d) difference image, frame (b) minus frame (a). Solar rotation 
compensated for all images to 00:40. 
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Figure 8. NoRH 17 GHz difference images (nonlinearly displayed). Black crosses mark the 
measured leading edge of the bright feature. White squares mark the measured center of the 
dark feature. Slanted cross marks the eruption center. 

at 17 GHz with an integration time of 1 min. Figure 8 shows some of them. 
Moving features arc still poorly visible, especially by cessation of the flare. Using 
extra averaging over three frames and subtraction of an image averaged during 
01:41-01:50, it is possible to detect two moving features in 17 GHz images. 

One of them is dark and moves toward the North Pole (white squares mark its 
center). Its lowest brightness temperature reaches 6 000-8 000 K against the quiet 
Sun's background (10 000 K at 17 GHz) pointing at its large optical thickness 
and, consequently, the kinetic temperature of ^ 6000 K. Another feature is bright 
and moves up in the plane of sky. Black crosses mark its roughly measured 
leading edge. It becomes pronounced after leaving the solar disk. This is probably 
due to its overlap with the dark feature, when both of them are not far from 
the eruption region. Thus, this feature is probably optically thin. Its brightness 
temperatures are 400-5000 K; hence, its kinetic temperature is probably » 
5 000 K. Its estimated velocity is ^ 500 kms^^, close to the speed of the bright 
feature Fl visible in TRACE 173 A images. Their trajectories are also close (see 
Figure 11). The bright NoRH feature resembles in shape the foremost bright 
envelope observed by TRACE at 00:15:20 (see Figure 3i). Although the NoRH 
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Figure 9. LASCO/C2 images of the related CME: ruiiiiiiig differences (a-d) and fixed-base 
differences {e-h, NW quadrant). All images are nonlinearly displayed. Thin circle: the solar 
limb, thick circle: occulting disk, small gray disk: the position of the eruption region, white 
dashed ovals outline the expanding CME. Labels denote discussed features. 



feature is not detectable during TRACE observations because of bright flare 
emission, these facts demonstrate the NoRH and TRACE bright features to be 
identical (some dissimilarities are due to their different temperature responses). 

2.6. CME in S0H0/LASC0/C2 Images 

LASC0/C2 & C3 observed a CME with a central position angle of 294° in 
many frames from 00:54 onwards. The CME was classified in the Preliminary 
2004 SOHO LASCO Coronal Mass Ejection List^ as a 'partial halo'. Figure 9 
presents four LASC0/C2 images. A bright part of the CME has an angular 
span < 90°. The spatial structure of the CME docs not match a classical three- 
component one. Its components do not resemble a self-similar structure with 
each one being inside another one. No cavity is detectable. The foremost part of 
the CME docs not resemble a well-defined frontal structure. It consists of diffuse, 
faint material, inside which an interconnected structure is visible (frame /) . 

The fastest spike-like feature 1 was used for measurements in the SOHO 
LASCO CME Catalog^; its linear-fit speed is 409 kms~^ Labels 2 and 3 mark 
the brightest components of the CME. Rotation of feature 3 determines its 
helical structure. The ejected material is detectable around almost the whole 
occultor, but it is faint in position angles — 270°, even relative to the foremost 
NW part (the fainter, faster SW part is listed in the SOHO LASCO CME Catalog 
as another CME). White dashed ovals in Figure 9 roughly outline the whole 
CME. Outflow of the ejected material continues a few hours {e.g., feature 4 in 
frame h). 



^ftp://Iasco6.nascom.nasa.gov/pub/lasco/status/LASCO_CME_List_2004 
^http:// cdaw.gsfc.nasa.gov/CME_list/ 
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Feature 3 appears from behind the occultor at a position angle corresponding 
to the trajectory of bright NoRH and TRACE features. Their shapes are rather 
similar. The position angle of feature 2 roughly corresponds to TRACE feature 
F2 and the bright sector visible in the EIT 195 A image at 00:24. 

2.7. Total Flux Radio Data 

Records of total radio flux made with Nobeyama Radio Polarimeters (NoRP) 
and in Learmonth and Ussuriysk Observatories show a strong burst to start at 
about 00:12:20 and reaching 1200 sfu at 17 GHz. Then, after 00:30-00:40, a post- 
burst decrease is observed at frequencies < 5 GHz by 01:15 01:50 (Figure 10). 
The duration of the 'negative burst' is minimal at 5 GHz and increases toward 
lower frequencies as well as the depression depth, from 5% at 5 GHz to 12% 
at 1 GHz. The deepest depression occurs at 00:55. The large-scale darkening is 
registered in the EIT 304 A channel at the final stage of the 'negative burst'. 

Thermal free-free emission of post-flare loops contributes to the total flux. The 
dotted lines in Figure 10 show the total flux records without this contribution. 
To estimate it, we computed emission measure and temperature, T, from soft X- 
ray GOES data following White, Thomas, and Schwartz (2005). The size of the 
soft X-ray emitting region was found from GOES/SXI images, and a geometrical 
depth, L ~ ^1/2 1.3- 10^ cm, estimated from its area, A. The optical thickness 
at a frequency u, = ki,L, was determined from an expression 

» 0.2nliy-'^T-^/'^L (CGS units). (1) 

2.8. Overall Picture of Near-Surface Activities 

Figure 11 presents an overall picture of large-scale disturbances observed on the 
solar disk and in its closest vicinity during and after the 13 July 2004 eruptive 
event. The figure summarizes observations in the Ha line, in 195 A, 304 A (EIT), 
and 173 A (TRACE) bands as weU as at 17 GHz (NoRH). 

The event starts at 00:02:30 with stretch of a filament or two 'merged' ones. 
During its steady rise, flare energy release starts by 00:13:40 beneath. The 
filament explosively erupts at about 00:15. The explosion produces Moreton 
and 'EIT' waves, which run faster than all other observed disturbances. Dark 
fragments of the filament form surges or fly mostly north, even reaching vicinities 
of the North Pole. One portion of the bright material moves up in the plane of 
sky, while another one moves west-northwest. The former portion is traceable in 
TRACE and NoRH images from the eruption site up to the edge of the field of 
view. The latter portion is only detectable in EIT 195 A images. 

3. Absorption 

As multi-spectral data show (Figure 11), fragments of the eruptive filament dis- 
persed over a large area could be responsible for absorption phenomena observed 
in radio range ('negative burst'), Ha (surges), and EUV at 195 A and 304 A. 
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Figure 10. Total flux time profiles in a range of 1—9.4 GHz recorded in Nobeyama, Learmonth, 
and Ussuriysk observatories. Dotted lines show subtracted free- free emission of the arcade 
computed from soft X-ray GOES data. 



In this section, we estimate parameters of the absorbing material and find out 
if it is possible to reconcile all mentioned manifestations. We then compare the 
estimated masses of the absorber and the CME. 

3.1. Simulation of Radio Absorption 

Wc simulate the 'negative burst' using a simple four-layer model with a uniform 
absorbing slab (Figure 12). The model contains 1) the chromoshpere, 2) the 
'spray' of area As 'inserted' into the corona, 3) a coronal layer of depth z 
between the chromosphere and spray, and 4) a coronal layer between the spray 
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Figure 11. An overall map of large-scale disturbances observed during the event on the solar 
disk and in its vicinity. Gray scale background: the EIT 195 A difference image (00:24-00:00). 
Blue contours: darkening in the EIT 304 A image (13 Jul, 01:20 - 12 Jul, 19:20); bright 
—20%, dark —60%. Green-yellow oval and arc: wave fronts in the EIT 195 A difference images 
(00:24 and 00:36). Black slanted cross: eruption center. Red crosses: measured positions of 
the Moreton wave fronts. Red arc with the arrow: approximate propagation direction of the 
wave along the solar surface; red arrows roughly show its angular span. Three brown arrows: 
trajectories of moving fragments, which were revealed from MLSO Ha images. Triangles: the 
leading edge of the expanding bright feature in the TRACE 173 A images. Circles: measured 
positions of the moving bright (white) and dark (black) features in NoRH 17 GHz images. 
Black-white dashed circle: the solar limb. 
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Figure 12. A simple model to estimate the frequency dependence of radio absorption. 



2004-07-13_preprint.tex; 6/11/2008; 9:18; p. 15 



(Tr('ciiii(*\" ct. al. 



and observer. The ratio of the observed flux to the quiet Sun's total flux is 

F/Fqs = [T<fs(A - As) + TiAs]/{T^sA). (2) 

A{i>) is the quiet Sun's area. We used frequency dependencies of the quiet Sun's 
radio radius and brightness temperature found by Borovik (1994). Superscripts 
'B' denote brightness temperatures to distinguish them from kinetic ones. Since 
the brightness temperature after each layer is a sum of its own emission and a 
non-absorbed remainder emission from preceding layers, the observed brightness 
temperature, T^, of the spray with a kinetic temperature Tg is 

Ti = Tchre-(^^+"'=+^^) + Tc(l - e-^i)e-(^=+^s) (3) 
+ Ts(l - e-^^)e-^= + Tc(l - e'^^), 

where Tj is determined by expression (1), ti = tc — T2, with rc being the 
total optical thickness of the corona at a given frequency, T2 = rce"^^/^, 
H = 2kTc/{migQ) k, 8.4 • 10^ cm the height of the uniform atmosphere, k 
Boltzmann constant, z the height of the absorbing layer above the chromosphere, 
rrii the average mass of ions, and Qq = GMq / Rq the gravity acceleration at the 
photosphere, G the gravity constant, Mq the mass of the Sun. 

Using this model, we fitted the radio absorption depth measured from multi- 
frequency records (Figure 10) for two instances, the maximum radio absorption 
(00:55) and EIT observation at 304 A (01:19). The results are shown in Fig- 
ure 13 with best-fit parameters. The model is rough; nevertheless, variations of 
parameters of ±20% cause detectable discrepancies with observations. 

With an optical thickness evaluated, a density and mass of the absorber could 
be estimated if its geometrical depth was known. The latter can be found from 
several facts. One estimate is the POS distance of the overturn of feature 3a, 
190 Mm (Section 2.1); with the position of the active region, it is expected to 
be comparable with its height. The time of the overturn, 00:38-00:40, provides a 
better estimate. Since the maximum height is large, we used general expressions 
with an inconstant g{R) = GMq/B? for a free fall time from the heliocentric 
distance r to R, velocity v at R, and r{v): 

* = \l 2g{R)R^ { + + a'^^^i^^l - 2^/^)] } ' 

v = y^2g{R)R{l-R/r), r ^ 2g(R)R^ /[v^ + 2g{R)R]. 

Using these expressions, we find the height of the overturn to be 180-200 Mm 
and the initial velocity of 350-380 kms~^, reasonably close to the measured 
POS parameters of the filament. From the height and time of the overturn, a 
height at 00:55 is 100-130 Mm. With z ~ 30 Mm estimated from the model, the 
geometrical depth is 70-100 Mm, and from expression (1) we find the density 
of (1.1 - 1.3) • 10** cm--"* and the mass of (3 - 4) • 10^'"^ g at that time. 

It is difficult to estimate the depth at 01:19 from observations. We speculate 
that the height at that time was significantly less than it at 00:55 was — e.g., 
< 50 Mm. On the other hand, there is no reason for the absorbing cloud to be a 
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Figure 13. Radio absorption measured directly (crosses) and with subtracted contribution 
from tiie arcade (circles) along with results of modeling (lines) for two instants. 



thin film located at a height of 10 Mm (from the model); therefore, we assume 
its depth to be > 2 Mm (most probable 10-20 Mm). With these depths, we get 
limits for the density of (0.5 — 2.3) • 10^ cm^"^ and mass of (0.8 — 3.5) • 10""^^ g; 
most probable are ~ 10® cm""^ and ~ 2 ■ 10^^ g. 

3.2. Absorption in Ha and EUV 

Ha opacity is due to neutral hydrogen absorbing continuum radiation of under- 
lying atmosphere. Mcin et al. (1996); Heinzel, Mein, and Mein (1999); Molowny- 
Horas et al. (1999), and Tziotziou et al. (2001) calculated physical parameters 
in stationary or slowly moving cool clouds like filaments on the basis of a 
non-LTE radiative transfer approach. Heinzel, Schmicdcr, and Tziotziou (2001) 
and Heinzel et al. (2003a) studied relations between opacities in Ha and Ly- 
man continuum at 912 A and proposed an approximate expression to estimate 
the hydrogen density from the optical thickness at the Ha line center. The 
Doppler-brightening effect decreases absorption of moving prominences, when 
line-of-sight velocities exceed ~ 50 kms~^. The Doppler shift of the Ha absorp- 
tion line exceeds the bandwidth of the BBSO Ha 0.5 A Lyot filter at a lower 
velocity of ~ 20 kms~^. 

A spectroscopic model for absorption of EUV emission in filaments and promi- 
nences was developed by Heinzel et al. (2003a, b), Anzer and Heinzel (2005), and 
Schwartz et al. (2006). They showed that the EUV opacity is due to photoion- 
ization of hydrogen (below the edge of the Lyman continuum, 912 A), neutral 
helium (< 504 A), and ionized helium (< 228 A). Absorption of continuum 
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Figure 14. The mass of the cjecta estimated from opacities in Ha (circles), 195 A (squares), 
304 A (triangle), and radio range (bars; limits, cross: probable). 



does not depend on radial or turbulent velocity. For 195 A, we take an average 
absorption cross-section of a 'cloud' consisting of 92% of H and 8% of He, 
fi95 « 7 • 10~^° cm^, which varies within 30-40% in a temperature range of 
< T < 80 000 K (we calculated the absorption cross sections for 195 and 
304 A using Tabic 2 and expressions (7) and (9) from Anzcr and Hcinzcl, 2005). 
Along with absorption, darkening in coronal lines could be due to a volume 
blocking effect (VBE), i.e., the absence of a coronal radiation from the volume 
of the cloud. This effect depends on the height and depth of the cloud relative 
to a height scale at a given wavelength. For 195 A, a typical radial height 
scale is ~ 70 Mm, and for our position angle, 42°, the line-of-sight height is 
H\\ ~ 100 Mm. 

We computed the hydrogen density from the optical thickness in the Ha 
line center assuming a line-of sight velocity V|| < 20 kms~^, a turbulent one 
Kurb = 5 kms~^, temperature 8 000 K, and depths of absorbing fragments to be 
equal to their widths. Figure 14 presents masses estimated from absorption of all 
emissions considered. The mass estimated from Ha increases from 00:14 by 00:30 
and then decreases likely because velocities of fragments become V\\ > 20 km s~^ . 

To estimate the opacity at 195 A, we analyzed EIT difference images (sub- 
trahend is 12 July, 23:59:57) in the interval of 00:12-02:00. Regions of probable 
absorption were found as depressions below la level, where cr was computed as 
an r.m.s. over negative pixels of darkening regions in the NW quadrant only. The 
regions revealed in this way contained long-lived dimmings {e.g., Dl, D2), which 
were unlikely to be due to absorption. We identified them from the difference 
image at 02:00 and excluded from computations. The mass estimated from 195 A 
images reaches a maximum of 3.7-10^^ g at 00:35 and then diminishes. The esti- 
mations miss some fragments of the cjecta that crossed the long-lived dimmings 
and could contain 20-30% of the total mass. A noticeable part of the ejected 
mass could be contained in hot ionized plasmas of the transition layer between 
cold absorbing material and ambient corona (analogous to an EUV filament 
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channel) seen in TRACE 173 A as a bright contour. We did not consider the 
VBE, because the depth of absorbing fragments was ^ H\\ . 

A large depression area at 304 A (01:19, sec Figure 7) shows a part of ejected 
material to return back onto the Sun. The same selection criteria as for 195 A 
were applied. A mean brightness ratio to the pre-event image in the whole region 
is about 0.54, while in its significant part the ratio is noticeably below 0.5. 

Absorption at 304 A can be due to photoionization of H and Hei by continuum 
or due to resonant scattering of the Hen line in He ions. The latter mechanism 
cannot provide the absorption ratio < 0.5: the cloud emits resonant photons in 
47r angle, while incident photons come from below and lateral directions, i.e., 2-k. 
Also, the resonant scattering is very sensitive to the line-of-sight velocity: the 
Doppler shift exceeds a typical width of the Hell 304 A line if V|| > 50 kms~^. Ab- 
sorption of the continuum is more probable (cross section (7304 ~ 5.5-10"^® cm^), 
and the mass is ^ 2.8 • 10^''' g (triangle in Figure 14). 

The closest 171 A (01:00) and 195 A (01:25) images do not show absorption in 
the 304 A darkening region (absorption at 195 A occurs in other regions). At this 
time, the cloud was likely in between of the heights of Hen 304 A illumination 
and coronal loops responsible for the 171 and 195 A emissions, as absorbing 
flows in TRACE 173 A images show. Low clouds with a temperature close to 
the chromospheric one, 10 000 K, are not detectable in radio range also. 

The CME mass estimated using LASCO software reaches 1.3 • 10^^ g at 04:06. 
Correction due to the off-plane direction of the CME (Vourlidas et al., 2000) of 
~ 1.4 gives its real mass of 1.8 • 10-^^ g. A lesser mass of 1.8 • 10^^ g in the SOHO 
LASCO CME Catalog is related to the faint SW part (see Figure 9a~c). 

The estimates of the mass of the ejecta from absorption of different emissions 
agree with each other within their accuracy and limitations of methods. The 
masses of the CME and the returned part are close. Thus, the ejecta disintegrated 
into parts with comparable; niassc;s. one of which flew away, while anotlic;r one fell 
back onto the Sun. One might ask if the absorbing cloud left the Sun later and 
was observed as feature 4 in Figure 9h. However, its estimated mass is 1.1-10^* g, 
much less than the mass of the cloud. This rules out such a possibility. 

4. Kinematics 

In this section, we choose analytic expressions to describe kinematics of a wave 
expanding from the eruption site and mass fragments flying away. Their birth in 
a single process is certain, but their subsequent interrelation is vague, i.e., it is 
not clear if a piston existed to determine the motion of the wave. The oval shape 
of the wave front disfavors the role of the faint NW sector in the EIT image 
at 00:24 as a piston. The wave and fragments appear to move independently. 
Deceleration of mass fragments, constituting the CME, points at importance of 
gravity. This also shows up in the fact that a part of the ejected mass falls back 
and produces a cool absorbing 'screen' to cover a significant part of the Sun. 

4.1. Approach 

We use a self-similar approximation of the expansion of mass fragments. Low 
(1982, 1984a,b) first applied it to expansion of a CME. Uralov, Grechnev, and 
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Figure 15. Distance-time plots (in solar radii) of fragments observed with TRACE at 173 A 
(gray erosses), NoRH at 17 GHz (circles), and LASCO/C2 (black; slanted crosses for feature 
1, straight crosses for feature 2, and triangles for feature 3). The large circle denotes the 
observation point of the bright aureole in the EIT 195 A image at 00;24. Dash-doted line 
denotes the position of the eruption center. 



Hudson (2005) derived expressions convenient to compare with observations that 
apply independent of the presence of a CME-driven coronal shock. If it is present, 
then CME parts are assumed to be confined by a contact surface, which is also a 
surface of a moving piston. A solution belongs to a class of self-similar motions, 
in which Lagrangian velocities of pieces of mass v are proportional to their 
distances r from the center of expansion, w oc r (see also Sedov, 1981, p. 318). 

To describe wave-like phenomena, we employ concepts of self-similar theory 
of a point-like explosion (e.^., Sedov, 1981; Zel'dovich and Raizer, 1966). Such 
an explosion produces a strong shock wave in the ambient medium. 

A crude example to combine both these approaches is an air explosion of 
a bomb consisting from a charge and a shell. The outcome is a shock wave 
and fragments of the shell flying away. They expand synchronously during the 
explosion and move independently afterwards. 

4.2. Mass Fragments 

As Section 2.8 shows, some counterparts of CME features are detectable on 
the Sun or nearby. Figure 15 shows distance-time points measured in POS 
from TRACE, NoRH, and LASCO data. The best coverage has feature 3 (see 
Section 2.6) and its counterparts. A quasi-periodical component in plots of 
all LASCO features is due to their rotation. It complicates measurements and 
decreases their accuracy. 

Expression (A8), = Vq -I- 2^/ro (1 — ro/r), from Uralov, Grechnev, and 
Hudson (2005) with r ^ oo, w ^ Woo and vq = v{r = Tq) gives 2A/rQ = v"^ — v^: 

v-' = vl + {vl,-vl){l-ro/r). (5) 
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This expression shows that with the increasing distance, the deceleration of all 
features monotonically decreases up to zero. Therefore, no polynomial fit of the 
measured data points applies. To compare results with observational data, a few 
parameters are necessary; the initial and asymptotical velocities, vq and Voo] the 
initial size of the ejecta, ro; and origins of measurements both in space and time. 
The origins are known with a satisfactory accuracy (see Section 4.4; start time of 
00:14:50). The t;oo can be estimated from Figure 15 (it is 409 kms~^ for feature 
1 according to the SOHO LASCO CME Catalog). The uo for feature 3 can be 
found from TRACE and NoRH data points. Obviously, w^^' > > 4'^'- The 
ro parameters could be only found in attempts to obtain a good fit. The curves 
in Figure 15 were obtained with the following parameters: 

LASCO feature 1: ro = 270 Mm, vq = 1100 kms-\ v^o = 409 kms'^; 

LASCO feature 2: ro = 220 Mm, vq = 850 kms-\ Woo = 310 kms'^; 

LASCO feature 3: ro = 110 Mm, va = 490 kms^^, = 100 kms'^. 

These parameters provided a satisfactory fit for each of the three features 
visible in LASCO images, with the accuracy being worse for features 1 and 2, 
whose counterparts are poorly identified at/near the Sun. 

4.3. Wave - Uchida's Model: Weak Fast Mode 

The almost perfect oval shape of the brightening in the EIT 195 A image at 00:24 
and its expansion visible in the next image at 00:36 suggest its wave-like origin 
and association with the Moreton wave. The following facts also support this: 

(i) its configuration is different from the main CME components, in particular, 

(ii) its eastern part is roughly as bright as the western part, if not exceeds it, 
whereas the eastern CME part in LASC0/C2 images is much fainter than its 
western part, and (iii) the apparent velocity of the oval, ~ 700 kms"-"^, is com- 
parable with the Moreton wave's one. This possibility was previously discussed 
by, e.g., Hudson and Warmuth (2004) and Warmuth et al. (2001, 2004a, 2004b). 

According to the original Uchida's (1968) idea, Moreton waves are excited by 
coronal waves, being their skirts. Uchida (1968) considered a weak short fast- 
mode MHD wave to propagate in a radially- symmetric corona, with the Alfven 
speed increasing with height. This caused refraction of rays into regions of lower 
Alfven velocity, toward the solar surface, and turned the wave front down. This 
model predicted acceleration of Moreton waves. However, Warmuth et al. (2001, 
2004a) found their systematic deceleration (Yamaguchi et al. (2003) reported 
acceleration of a Moreton wave in its propagation into a coronal hole). Their 
observed velocities are sometimes too high for weak fast-mode MHD waves. 
These inconsistencies with observations do not appear to rule out weak fast- 
mode MHD waves completely, but encourage to search for another possibility 
to match observations better. Warmuth et al. (2001, 2004a, 2004b) and Hudson 
and Warmuth (2004) proposed blast waves to be responsible for Moreton waves 
and, probably, some 'EIT waves'. 

4.4. Wave - Strong Shock 

We use estimative expressions from a theory of a strong point-like explosion in 
a variable-density medium. A rigorous analysis of this gas dynamic task was 
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carried out by Sedov (1981). Self-similarity of the solution is ensured by a large 
pressure excess inside the volume confined by the shock front over non-disturbed 
medium. If the magnetic field in ambient plasma is strong, then the profile and 
shape of the shock front change. This change is small for a strong shock, because 
the increase of the magnetic pressure in the shock is determined by extreme 
plasma compression, whereas gas pressure is unlimited. Also, the magnetic field 
in our case is weak, because the disturbance propagates over quiet Sun's regions. 

Kinematics of a single-pulse spherical blast wave is governed mainly by pile- 
up mass. We consider propagation of a self-similar blast wave excited by an 
explosion of an energy E in media (a) with constant density, and (b) with a 
radial density falloff from the explosion center, p oc r~" . 

1. For a constant density, p = po'. poR^V^ = const E [R radius and V velocity 
of the shock front), with const corresponding to constancy of potential to 
kinetic energy ratio for plasma involved into self-similar motion. Thus, V = 

^ = (2^^^) ^'^ oc and R oc e'K 

2. Similarly, for p = br'"': f ■^r'^dr l/^^P-a] = ^^^g^ ^ Consequently, 
F=^oci?-K3-«)/21,and 

i?octP/(5-«)]=t5. (6) 

A strong spherical shock decelerates if a < 3 and accelerates if a > 3. Thus, 
5 = 2/(5 — a), with S = 0.4 for a homogeneous medium {a = 0) and S = 
0.67 for a medium with p oc r~^. To compare expression (6) with measured 
points of the Moreton wave, we use expression (6) in a form {R — Rq)/{R — 
R\) = [{t — to)/(t — ti)]'^^^^~"-". Here Rq is the distance between the origin of 
measurements and explosion center, to the corresponding time, and ti, R\ are 
related to one of measured points. 

Figure 16 shows the results. Triangles denote the Moreton wave, straight 
crosses denote distances of the EIT oval from the eruption center measured in 
the same direction, slanted crosses show measurements of the EIT oval toward 
the solar disk center, and a rectangle presents the first manifestation of the wave 
m a TRACE 173 A image (width: exposure time, height: positional uncertainty). 
The Moreton wave is expected to slightly lag behind the coronal wave (Vrsnak 
et al, 2002), as the TRACE and Ha fronts show. We used the last point of the 
Moreton wave as a reference one (ti,Ri). The dashed curve corresponding to 
a = shows much stronger deceleration than actually observed. A good fit is 
achieved with a — 2 {5 = 0.67) and a start time of 00:14:50±20 s for both the 
Moreton and 'EIT' waves. The positional uncertainty of the explosion center is 
< 20 Mm; origins of both the position and time are close to the observed eruption 
ones. The difference between distances measured along the are on the sphere in 
Figure 11 and POS distances in the same direction does not exceed measurement 
errors. 

The kinematical closeness of the Moreton and 'EIT' waves confirms their 
common nature, and agreement with the calculated plot supports their origin 
due to coronal blast shock. Note that data points of the 'EIT wave' suggest a 
stronger deceleration, implying a density falloff a < 2 at longer distances. 
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Figure 16. (a) POS Distance-time plots of the Moreton wave (triangles) and EIT oval 
(crosses) measured in the same direction from the eruption center. Rectangle denotes the 
first manifestation of a wave in a TRACE 173 A image. Curves: strong shock with density 
distributions of p = const (dashed) and p oc r^^ (solid); weak shock with p = const (dotted). 
Slanted crosses: distances of the EIT oval from the eruption center measured toward the solar 
disk center; dash-dotted curve shows their fit with p oc r~^. (b) Instantaneous velocities 
computed for a strong shock with p oc (solid) and a weak shock with p = const (dotted). 

The calculated instant shock speed in Figure 16b can be compared with 
estimates from type II bursts listed in Solar Geophysical Data"^. The estimates 
decrease from 1000 kms-^ at 00:17-00:20 (Culgoora) via 850 kms"^ (00:16- 
00:40, Learmonth) to 550 kms~^ at 00:27-00:41 (Holloman). Such estimations 
use model plasma densities, which are uncertain. Nevertheless, the calculated 
shock speed and estimations from type II bursts are in overall agreement. This 
supports the same origin of the Moreton wave, 'EIT wave', and the type II 
radio burst, likely associated with the same decelerating coronal blast shock, as 
Warmuth et al. (2001, 2004a, 2004b) concluded (however, Pohjolainen, Hori, and 
Sakurai (2008) interpret this type II burst to be due to two distinct shocks). 

Our plot in Figure 16 was obtained under assumptions of (1) a strong shock 
and (2) an omnidirectional density falloff from the eruption site. Their agreement 
with observations appears to be surprising, because a horizontal density falloff 
is not expectable at large distances. We discuss this issue in the next section. 

We also fitted the oval envelope of the CME with p oc r"^ '' expectable for 
a coronal shock propagating outward: at heights of (0.2 — lO)i?0, the density 
in the corona above an active region falls off oc (r the distance from the 
photosphere) according to Newkirk model (Newkirk, 1961), and p oc r~^'^ from 
a model compiled by Gary (2001). This fit is close to the motion of LASCO 
feature 1. Note that with <gC Vq and a = 2, expressions (5) and (6) predict the 



^http://sgd.ngdc.noaa.gov/sgd/jsp/solaxindex.jsp 
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same height-time plots. Poor measurement accuracy does not allow to distinguish 
between the wave and mass fragments. Their closeness is also possible. 

4.5. Comments on a Realistic Situation 

We come to the following possibilities for the Morcton wave. 

1. Uchida's (1968) refraction model is a linear acoustic approximation, which 
is very sensitive to variations of plasma parameters. This model predicts accel- 
eration of a wave that docs not correspond to the observed situation. 

2. Strong shock wave: gas velocity behind the wave front (Lagrangian velocity; 
it is presented by subsequent motion of loop LI m TRACE images, Figure 3, 
350 kms"'^) is much higher than the ambient fast magnctosound speed, Vf. This 
strongly non-linear wave is mainly sensitive to density distribution. Deceleration 
naturally appears and critically differs if the wave runs along, against, or across 
the density gradient. To explain propagation of the Morcton wave, we assumed 
a spherically-symmetric coronal density falloff from the explosion center. In this 
case, the wave runs contrarily to the density gradient, and its deceleration is 
determined by expression (6). This resulted in a good fit with a = 2. Otherwise, 
the deceleration is too strong, as the following considerations show. 

The Moreton wave runs horizontally, presumably being a skirt of a coronal 
shock front. With a constant density at a fixed height, a strong shock decelerates 
even faster than expression (6) with a = predicts. To roughly estimate this 
case, we consider a vertically-stratified atmosphere with an upward density falloff 
of p oc r~". Here, the coronal part of the Moreton wave runs across the density 
gradient. The front of a blast wave in such atmosphere is not a spherical one; 
its initially vertical parts in propagation progressively incline down at the same 
height. We assume pressure P to be uniform inside a volume confined by the 
shock front, P oc with R being the vertical extension of the wave surface. Its 
variation is determined by expression (6). In a strong shock, PxV^ ~ P, where 
Vx = '^^j^ is a velocity of the front along a constant-density surface, px = const. 
With a = 2,Vx oc t~^ , and Rx oc ln(t). After the fast onset, such a shock rapidly 
decelerates up to Vf (cf. a typical speed of 'EIT waves'). 

Besides a spherically-symmetric density falloff, a high horizontal deceleration 
of a strong shock could be prevented if a piston moved properly in the same 
direction (cf. Low, 1984a). If a self-similar solution does exist, then a piston and 
wave of a common origin diverge slowly. Indeed, the velocity of the strong shock 
front is ^ = ^^^Vgh with 7 the polytropic index and Vsh the Lagrangian plasma 
velocity just behind the shock front. The velocity of the contact discontinuity 
is ^ = Vc with Vc the Lagrangian plasma velocity within the discontinuity. In 
such self-similar motions, ^ = From these expressions, ^ = and 
the ratio of distances of the shock front R and the piston from the explosion 
center is ^ oc ri'^~^^^'^. Initially R = rc = I, and then the ratio increases slowly; 
both the piston and wave are expectable to be not far from each other in images, 
which we showed. However, the wave appears to have run far away well before 
the first appearance of mass behind it. A blast wave is also favored because the 
initial velocities of the ejecta do not point in the same direction as the EIT oval. 
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An omnidirectional density fallofF seems to be possible near the active region, 
but unlikely at long distances. A possible solution of the problem is subsequent 
damping of a strong shock to a moderate or weak one. The weak shock runs 
along rays determined by the spatial distribution of Vf, which, in turn, depends 
on the density distribution. The speed of the weak shock depends on the density 
gradient along the ray {e.g., horizontal) rather than on the crosswise (vertical) 
gradient. Just the vertical density gradient caused the extreme deceleration of a 
strong shock propagating horizontally, that was considered above. 

3. Weak shock: gas velocity behind the wave front is < Vf. Propagation of a 
spherical weak shock in a uniform plasma calculated using expressions from a 
paper of Uralova and Uralov (1994) is shown in Figure 16 by dotted lines for 
V{ = 300 kms"'^ and the phase speed of the wave at its first manifestation of 
1600 kms~^. This curve matches the coronal wave, whereas at shorter distances 
its agreement with the Moreton wave is worse. The initial speed is too high for 
the approximation of a weak shock, suggesting a medium-intensity shock. 

4. Shock of a moderate intensity. This intermediate case is difficult to calcu- 
late. The curves in Figure 16 imply that the wave is initially a rather strong one, 
and, after some transitional stage, it becomes weak. 

5. Large- amplitude simple wave (without discontinuity). We calculated this 
case for a uniform density distribution formally using a solution for a simple 
wave and a spherical wave front. The simple wave has a distance-time plot 
close to the weak shock fit in Figure 16; its velocity is less at short distances. 
However, the solution of a simple wave rapidly becomes incorrect. Estimations 
show that the shock discontinuity appears either during formation stage of the 
wave disturbance, or just afterwards. To fit the calculated distance-time plot to 
experimental data points in Figure 16, the maximum Lagrangian velocity must 
be a few times higher than Vf. The appearance of a type II burst as early as at 
00:16-00:17 was probably due to the rapid steepening the wave into a shock. 

All these curves match more or less observations. Data available do not to 
permit us to find out a particular type of a wave. The spherical strong shock in 
decreasing density matches observations better and seems to apply at moderate 
distances. We therefore use it in our subsequent discussion of the Moreton wave. 

The fact that the best fit for the Moreton and 'EIT' waves corresponds to 
a = 2 is probably accidental. Our idealized assumptions certainly affect results. 
We assumed a strong blast wave and a radial density falloff from the explosion 
center, which in reality might be anisotropic and probably variable. The shape 
of the wave front is also important; the pile- up mass grows with distance in a 
spherical wave faster than in a cylindrical or a flat wave. 

Some range of parameters is expectable in observations. Warmuth et al. (2001, 
2004a, 2004b) found average power-law indices 5 for Moreton waves' distance- 
time plots of 0.57-0.62 corresponding to density falloffs of a = 1.49 — 1.77 for 
spherical waves. The whole range of the measured indices was mainly from 0.41 
to 0.91 {a = 0.12 - 2.8) with two exceptions of 0.29 and 0.34. Both last events 
were located on the limb that made correct identification of wave fronts difficult; 
however, such strong damping and deceleration were, in principle, possible if 
explosions occurred in media with uniform horizontal density distributions. 
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The power- law distance-time relations found by Warmuth et al. (2001, 2004a, 
2004b) and in our event are expectable for coronal blast waves, in accord with a 
result of Balasubramaniam, Pevtsov, and Neidig (2007) and the original idea of 
Uchida (1968) in a general form. Excitation of a filament by a Moreton wave from 
its upper edge in the 6 December 2006 event also could be due to an incident 
coronal front inclined down (Gilbert et al., 2007). The nature of Moreton waves 
due to self-similar coronal shocks is supported by linear growths of wave front 
thicknesses with distances found by Warmuth et al. (2001, 2004b, 2005) in some 
events. It is also supported by correspondence with LASCO data. 

5. Discussion and Conclusion 
5.1. Overall Scenario of the Event 

The event started with a steady rise of a filament in the active region. Then 
an explosive eruption occurred probably due to mergence of two close filaments. 
Uralov et al. (2002) argued explosive development of MHD instability in such 
situation; Hansen, Tripathi, and Bellan (2004) confirmed this experimentally. 
The eruption was likely driven by MHD forces, because the filaments remained 
cool. The fastest observed feature with a POS acceleration > 8 kms~^ was 
bright, which means that, at least, its parts were at coronal temperatures. 

The initial energy of the cjccta was partially spent to disperse the filaments 
over a large area. The remainder energy was insufficient for all fragments to 
overcome gravity; both escaping and returning features were observed. The im- 
portance of gravity is manifest in deceleration of all observed features. Since the 
ejecta disintegrated, the CME did not resemble a classical three-part structure 
with a self-similar frontal structure, cavity, and core. One part of the ejecta 
escaped as the CME, and another one fell back onto the Sun. Masses of both parts 
were close. The latter part consisted of fragments of filament (s) dispersed into a 
cloud covering almost the whole NW quadrant. The cloud absorbed background 
emissions that was observed as moving EUV dimmings and a 'negative radio 
burst'. The cloud later fell back onto the solar surface, slipping along magnetic 
field lines. Other fragments produced surges visible in Ha up to heights of ~ 200 
Mm. The initial speeds of the fastest parts of filaments were less than the speed of 
their bright envelope (~ 500 kms~^ in the plane of sky) by 30-50%; fragments 
of their legs moved still slower. The cooler material after the explosion was 
affected by gravity only, had not exceeded the escape velocity (618 kms^"'^), and 
mainly fell back. All observed returning features were cool, while no falling hotter 
material was seen. By contrast, all observed escaping features were significantly 
hotter, at about coronal temperatures, and just they had coimterparts in LASCO 
images. Since masses of escaping and falling parts were comparable, the CME 
was mainly at ~ 1 MK, while the returned part was mostly cool, at ~ 10* K. 

The explosion produced a decelerating coronal blast shock. The presence of 
a wave is demonstrated by a displacement of loop LI in Figure 3i,j, when the 
ejecta was far from it. The skirt of the shock was seen as a Moreton wave 
decelerating from the initial speed of > 1000 to ~ 600 kms~^, while it was 
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observed. At longer distances, the shock was observed as an 'EIT wave' with a 
speed of < 500 kms~^. 

The initial magnetic structure was closed and confined inside the filament (s). 
After the eruption, dispersed fragments flew along lengthy field lines toward the 
North Pole and landed on a huge area. Some magnetic structure was carried 
away by the CME, as rotation of its component demonstrates. It is difficult to 
reconcile these facts without magnetic reconnection, but its role in the event 
is not quite clear. Probably, just the destruction of the magnetic configuration 
could cause dispersing the cjccta rather than its ejection as a whole. Anyway, 
reconnection does not seem to be directly related with the observed 'EIT wave' 
and dimmings. 

5.2. Comments on Nature of Large-Scale Disturbances 

The bright oval 'EIT wave' in this event was likely due to the coronal blast shock 
as follows from its kinematics, which agrees with the motion of the Moreton wave 
and theoretically expected propagation of the shock. The spheroidal shape of the 
wave front and correspondence to parameters of the type II radio burst confirm 
this conclusion. Quantitative parameters inferred in our analysis and results of 
other authors obtained for several 'Moreton & EIT wave events' show them to be 
mainly consistent with their origin due to coronal blast shocks. They are proba- 
bly moderate-intensity ones or strong shocks subsequently damping to moderate 
waves. Shocks compress (hence brighten) and heat plasmas. Manifestations of 
heating were presumably found in observations of some 'EIT waves' (Grechnev 
et al, 2005; see also Warmuth, Mann, and Aurass, 2005). 

Other large-scale transients of low brightness (visible in Figure 5b north, 
north-northeast, east-southeast of the eruption center, and a sector above the 
NW limb) moved slower, had different shapes, and were probably due to a 
different reason. On-disk brightenings could be similar to the off-limb feature. 
However, they could be due to some processes in low coronal layers, i.e., different 
from the off-limb feature. Their nature is not clear. 

Two kinds of dimming were observed. Shallow dimmings moving behind the 
oval 'EIT wave' front were due to absorption of the Sun's background emission 
in the moving 'cloud' of filament fragments, as quantitatively confirmed. By 
contrast, deep, long-lived dimmings Dl, D2 in Figures 4,6 were due to significant 
changes of previously bright coronal loops. Figure 4 shows that dimming D2 was 
too large to be due to displacement of loops. A plasma density decrease is the 
most probable (if not the only) reason that can occur in two cases, (a) Previously 
closed loops open or stretch into interplanetary space that results in plasma 
outfiows. (b) Loops slightly stretch, but do not significantly change. As can be 
seen, with the same number of particles in a loop, a change of its volume from 
Vb to Vi correspondingly changes emission measure, EMi/ EMq = Vq/Vi. Thus, 
stretch of a loop decreases its brightness. In summary, deep, long-lived dimmings 
are likely due to density decrease in coronal structures and obviously cannot 
travel. These conclusions agree with those of Chertok and Grechnev (2005a). 
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5.3. Estimations of Mass from Absorption 

We estimated the mass of the ejected 'cloud' from absorption in the radio range, 
in the Ha hne, and in EUV lines of 195 A and 304 A. All these methods supplied 
close results, from which we conclude that the mass of the 'cloud' was ~ 3-10^^ g, 
a typical mass of a filament. The mass of the CME was ~ 2 • 10"'^^ g, also typical 
of CMEs. This implies that such a spectacular event was unlikely an exceptional 
one. On the other hand, consistency of the estimates obtained by means of these 
methods demonstrates that they offer a promising way to evaluate masses of 
on-disk ejecta. 

Each of these methods has advantages and limitations. Ha images distinctly 
show absorbing features, but estimations are constrained to low temperatures 
and velocities and need unknown geometrical depths. Estimations from 195 A 
emission are not velocity dependent, but they do not say much about the temper- 
ature of the absorber and fail at low heights. Estimations from the 304 A emission 
are sensitive up to the transition region, but crucially depend on temperature 
and absorption mechanism (velocity dependence is important for the resonant 
scattering). Estimations from multi- frequency radio data need unknown geomet- 
rical depth and fail at low heights, but they allow to estimate the temperature 
of ejecta in a wide range and do not depend on velocity. 

5.4. Summary and Concluding Remarks 

From multi-spectral data, an overall picture of the 13 July 2004 eruptive event 
was reconstructed and confirmed by quantitative estimations. An explosive fil- 
ament eruption occurred in an active region. The ejecta disintegrated into two 
parts of comparable masses, one of which flew away as a decelerating CME, and 
another part dispersed over almost the whole NW quadrant of the visible solar 
disk. The latter part absorbed background solar emission that was observed as 
widespread faint moving dimmings at 195 A, a 'negative' radio burst, and a huge 
dimming at 304 A. By contrast, deep, quasi-stationary dimmings also observed 
at 195 A in this event were likely due to density decrease in coronal structures. 

Properties of both the Moreton wave and the oval faint front observed at 
195 A in this event are consistent with their origin due to a coronal blast shock 
produced by the eruption. Our results show what could be reasons for 'EIT 
waves' and dimmings in a particular event, not pretending to apply to all events. 

The coronal blast shock, the Moreton wave, and a huge cloud of the dispersed 
material of filaments were interconnected in our event. It was remarkable, but 
probably not an exceptional one. CMEs produced in such events are not expected 
to accelerate. A polynomial fit is inappropriate to analyze their kinematics. 

Important information about eruptive phenomena could be obtained using 
the following methods. Radio observations in a range of 1-10 GHz can supply 
parameters of ejecta (multi-frequency imaging observations would be especially 
valuable). Space-borne 304 A and Ha observations with an interval of ~ 1 min 
would be important. Masses of on-disk ejecta can be estimated from absorption. 
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